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ABSTRACT 

Context. Millisecond pulsars are very useful for determining the properties of neutron stars, for testing General Relativity and for 
detecting gravitational waves. However, the number of known millisecond pulsars is very small compared to that of ordinary pulsars. 
Aims. We conducted one blind and three targeted searches for millisecond and submillisecond pulsars at radio frequencies. 
Methods. The blind search was conducted within 3° of the Galactic plane and at longitudes between 20° and 110°. It takes 22073 
pointings to cover this region, and 5487 different positions in the sky (i.e. 25% of the total) were actually observed. The first targeted 
search was aimed at Galactic globular clusters, the second one at 24 bright polarized and pointlike radiosources with steep spectra, 
and the third at 65 faint polarized and pointlike radiosources. The observations were conducted at the large radiotelescope of Nancay 
Observatory, at a frequency near 1400 MHz, the exact value depending on the backend. Two successive backends were used, first a 
VLBI S2 system, second a digital acquisition board and a PC with large storage capacity sampling the signal at 50 Mb/s on one bit, 
over a 24-MHz band and in one polarization. The bandwidth of acquisition of the second backend was later increased to 48 MHz 
and the sampling rate to 100 Mb/s. The survey used the three successive setups, with respective sensitivities of 3.5, 2.2, and 1.7 mjy. 
The targeted- search data were obtained with the third setup and reduced with a method based on the Hough transform, yielding a 
sensitivity of 0.9 mJy. The processing of the data was done in slightly differed time by soft-correlation in all cases. 
Results. No new short-period millisecond pulsars were discovered in the different searches. To better understand the null result of 
the blind survey, we estimate the probability of detecting one or more short-period pulsars among a given Galactic population of 
synthetic pulsars with our setup: 25% for the actual incomplete survey and 79% if we had completed the whole survey with a uniform 
nominal sensitivity of 1.7 mjy. The alternative of surveying a smaller, presumably more densely populated, region with a higher 
sensitivity would have a low return and would be impractical at a transit instrument. 

Conclusions. The null result of our blind survey is compatible with our present understanding of the Galactic population of 
millisecond pulsars. In particular, there does not seem to be a large population of fast-rotating millisecond pulsars. 
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OO 1. Introduction tributed sample of millisecond pulsars has the potential to de- 

, tect gravitational waves at nanohertz frequencies (van Haasteren 

, In the past fifteen years the timing of millisecond pulsars in et al 2009 Verbiest et al Manchester!!^- Finally, the 

r-H , the radio domain has become an important field of research for discovery of subm illisecond pulsars would revolutionize funda- 

> ; pulsars and their emission mechanisms. Important results have mental phys[ ^ becaufje the stiffefjt equations of state would then 

^ , been obtained in neutron-star physics and General Relativity haye tQ be abandoned (H aensel et al.[T995l 
>^ ■ (see Lattimer & Prakash 2007 and Lorimer 2008 for reviews). 

U 1 It has also been shown that, for certain pulsars, the stability of 11 is likel y that onl y a sma11 fraction of the Galactic pop- 

the rotation period is close to that of atomic time scales for an av- ulatlon of millisecond pulsars has so far been discovered. One 

eraging time of several years (e.g. Janssen et al. l20T0l) . However, reason for a low detection rate is that the sensitivity of many sur- 

this property generally does not hold when the observations ex- ve y s decreases strongly with decreasing rotation period, mainly 

tend over decades. Nevertheless, assuming a good stability of the for hi S h values of interstellar dispersion, because of insuffi- 

rotation, a "pulsar" time scale may be calculated from a pulsar Clent sampling resolution in time and frequency. In particular, 

timing array (Tavlor [T99T1 Petit & Tavella [T996l Guinot & Arias the searches at the turn of the century for submillisecond_pul- 

I2QQ5I Manchester l201U sars were unsuccessful (D'Amico 120001 Crawford et al. 120001 

It is thus essential to discover new millisecond pulsars in Edwards et al.| 2001| l. 
order to improve on the present results. In addition to better This paper reports on a blind survey in the Galactic plane and 

statistics on a larger set of data, one can hope to discover pul- three successive targeted searches conducted at the large radiote- 

sars with a better stability of rotation. Presently, it seems that lescope of Nangay Observatory for short-period millisecond pul- 

old millisecond pulsars with short period (less than 15 ms) and sars (ranging from 0.65 to 30 ms). The project was initiated in 

thus small period derivative are the most plausible candidates the 1990s to search for short-period pulsars with a technique 

for a stable period. Furthermore, the timing of a widely dis- of soft-correlation of radio signals developed for other applica- 
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tions. At the time, the recording of radio signals at a very high 
sampling rate and their subsequent treatment by soft-correlation 
techniques was quite a challenge in view of the performances of 
standard computers, and was only being done in VLBI with ded- 
icated instruments. The first years were devoted to writing soft- 
ware to read the VHS (and later SVHS) video tapes on which the 
data were initially stored and to developing dedicated parallel 
computers based on transputers, and later on shark processors, 
to perform the actual soft-correlations, along with the necessary 
software. The blind survey started in 1997, a time when past 
surveys had searched only a fraction of the parameter space, be- 
cause of the limited sampling resolution in time and frequency. It 
was designed to be complementary to past and ongoing surveys 
by only searching with high sampling resolution. When the blind 
survey proved less successful than expected, targeted searches 
were meant to provide additional chances of finding pulsars with 
limited additional observing time. 

The acquisition method is explained in Sect. [2] the data pro- 
cessing in Sect. [3] and its validation in Sect. [4] The blind survey 
and the three targeted searches are presented in Sect. [5] We put 
the null result of the blind survey in perspective by estimating 
in Sect. |6]the chances of detecting short-period pulsars with our 
setup, by way of simulations on a model for the Galactic popu- 
lation of millisecond pulsars. 

2. Method of acquisition 

The observations were obtained at the large radiotelescop43 of 
Nancay Observatory. This is a transit instrument with a beam of 
3.6x22 arcmin at zero declination and 1400 MHz. 

In the initial stages of the project (1997-98), the data acquisi- 
tion was done with the VLBI S2 system and SVHS video tapes. 
The signal from the receiver in two bands (1400 and 1410 MHz) 
and in both polarizations was registered, totalling four channels, 
each 6.4 MHz wide. The sensitivity of the system was deter- 
mined from the mean flux density of a detectable pulsar (7 sig- 
mas). We measured a sensitivity of 3.5 mJy for 240 s of inte- 
gration, from the signal -to-noise ratio (SNR) of PSR B1937+21 
(Rougeaux et al. 2000). We obtained 186 positions in the sky 
(hereafter points) of the blind survey with this setup. Shortly 
thereafter, the radiotelescope closed for renovation of the re- 
ceiver, and the French space agency discontinued its support for 
this project, thus bringing observations and technical develop- 
ments to a temporary halt. 

In 2004 we replaced this system with a custom-made acqui- 
sition board that digitized the signal on one bit before sending it 
to an electronic board ADLink PCT7300A revB, which stored 
the data on the 120Gb hard disk of a PC running under Windows 
2000 OS. The system sampled the signal at a rate of 50 Mb/s 
on one bit and a 24-MHz band, one polarization, one channel, at 
the frequency 1370 MHz. One channel of this system (24 MHz) 
was equivalent to the old one (4x6.4 MHz), but the new sys- 
tem was much simpler to operate. We also confirmed the gain in 
sensitivity by a factor 2.2 of the new receiver of the radiotele- 
scope after 2000. Procedures allowed the automatic acquisition 
of the data and their analysis in slightly differed time on PCs. 
The computation of the autocorrelation function with small time 
samples (typically 32 fj.s) was done on that computer before fur- 
ther treatment on more powerful computers (part of a cluster of 
98 PC boards, and later two fast dual-core computers operating 
under Linux OS). The integration time of a single pointing was 

1 http://www.obs-nancay.fr/en/index.php/instruments/radio- 
telescope 



set to 200 s and the corresponding sensitivity was 2.2 mJy. We 
obtained 2829 points of the blind survey with this setup. 

In July 2006, a new custom-made acquisition board was in- 
stalled at the radiotelescope. The signal was thereafter obtained 
at 1392 MHz, with a 48-MHz bandwidth and sampled at a rate 
of lOOMb/s, which raised the sensitivity by a factor 1.4. The in- 
tegration time of each pointing was reduced to 160 s, so the sen- 
sitivity was 1.7 mJy. All three targeted searches and 2542 points 
of the blind survey were done with this setup. 

The functioning of the backend was regularly checked by 
observations of PSR B 1937 +21, about twice a month. 



3. Data processing 

Our data processing scheme for 1-bit data has been described 
earlier (Rougeaux 1999 Rougeaux et al. 12000 ). Here we recall 
the algorithm used to detect pulsar signal in a single observation 
(Sect. 13. U . This processing method is hereafter called "survey 
mode". We then present the algorithm used to increase the sen- 
sitivity by using multiple consecutive observations (Sect. 13.21 1. 
which we call "cluster mode". The blind survey was reduced 
with the first method, while the reobservations of pulsar candi- 
dates in the blind survey and the three targeted searches were 
reduced with the other one. 



3.1. Processing in "survey mode" 

The first step of the processing is coherent folding of the signal. 
Then a search in a 3 -dimensional parameter space (period P, 
dispersion slope D, and acceleration parameter A) is performed. 
The range of trial periods was initially 0.65 - 3.0 ms. The lower 
limit corresponds to the fastest theoretically possible rotation pe- 
riod for a uniformly rotating neutron star (Haensel et al. 1 1995b . 
Following the recommendation of the radiotelescope time allo- 
cation committee, we extended this range to 0.65 - 30 ms when 
we started the targeted searches at the end of 2006, at a low com- 
putational cost (about 25%). We then also applied the extended 
range to the blind survey, thanks to the increased computational 
power available. 

For a given set of trial values P, D, A, the folding process 
yields an SNR value S (P, D, A). The typical numbers of trials are 
N P - 1 6 -2 xlO 6 for period, N D - 100-200 for dispersion slope, 
Na < 10 for acceleration. In addition, for each trial, the signal is 
searched for in about ten positions in phase so that the search 
space has about 10 10 elements. All SNR values higher than a 
given threshold (e.g. for an SNR threshold of 4.0, typically 10 6 
elements out of 10 10 ) are kept for further analysis. Because the 
acceleration parameters are relatively few, the results are exam- 
ined independently for each value of this parameter. The retained 
values Sa(P, D) are then grouped by peaks (clusters of adjacent 
points in the plane P - D). 

Once this basic treatment is done, all signals with maximum 
SNR higher than some threshold (typically 6.5) are extracted and 
examined. Except when radio-interferences are present, the pro- 
cessing typically yields one such candidate peak every three to 
five observations. For the candidate peaks, the raw data are rean- 
alyzed with slightly different trial values, centered on the candi- 
date values, and with finer steps in period, dispersion slope and 
acceleration. Candidates that maintain the SNR for varying trial 
values (about one per 50-100 observations) are reobserved. 
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3.2. Processing in "cluster mode" 



For No consecutive observations, the results are then No series 
S i(P,D,A), with i from 1 to No- For a pulsar in a binary sys- 
tem, P may vary within each observation (i.e. the acceleration 
parameter is not zero), and P will also vary between each obser- 
vation, therefore a new parameter Pd - P(i + 1) - P(i) is needed. 
Obviously the parameters A, and Pd are correlated. However, 
the Aj are poorly determined in each observation, thus they can- 
not directly be used to connect candidate peaks between the No 
consecutive observations, especially as the candidate peaks are 
quite numerous (all SNR greater than typically 4-4.5). On the 
other hand, the dispersion measure of a pulsar is the same in all 
observations. The processing is then broken into No indepen- 
dent procedures, each concerning one trial dispersion slope, for 
which the matching of (P,A, Pj) is performed as follows. 



80sec, 24MHz-band 



The Hough transformation (Aulbert 2007) is used to trans- 
form the search in the plane (T,P), where T represents time 
(i.e. the sequence of observations), into one in the plane (P4, P). 
Assuming that the total time span is short enough for a pulsar 
period to vary linearly with time, in the first parameter space, a 
pulsar is represented by a line P = POd-i + P0 where PQ is the 
period at the reference epoch, and POd the period derivative. The 
line representing a pulsar is to be found as one that fits the mea- 
surement results (i, Pi). In the second parameter space (Pd, P), 
the measurement results are transformed to lines P = Pj - i.Pd 
and the pulsar is represented by a single point (POd, P0). The 
point representing the pulsar is found at an intersection of lines. 

With respect to this basic procedure (see a more detailed 
description in Aulbert, 2007), two additional features are intro- 
duced, as each measurement (i, Pi) comes with an SNR S , and 
an acceleration parameter A, . First, the acceleration parameter is 
used to constrain the P,/ value; i.e., each measurement produces 
not a full line, but a segment in the (P, Pd) plane. Second, each 
segment is assigned Si as a weight so that each intersection of 
lines comes with the combined weight of the participating seg- 
ments, in effect providing the detection SNR if the intersection 
of lines corresponds to a real pulsar. 

The results S {P0, POd) for all dispersion slopes D are then 
assembled and grouped by peaks (clusters of adjacent points in 
the plane P - D) and peaks with large SNR are examined. The 
couple (P0, POd) for the maximum SNR provides the period and 
period derivative of the observed, or candidate, pulsar. An exam- 
ple of such a processing is given in Sect.| 



4. Theoretical performance 

The point source efficiency at zero declination is about 1.6 K/Jy 
at 1400 MHz, and the system temperature is about 35 K, which 
leads to a 7-sigma theoretical detection threshold of 1.2 mJy 
for one 160-s observation at 100 Mb/s. The actual sensitivity 
is lower for several reasons: loss due to the 1-bit digitization, de- 
crease in sensitivity with declination differing from 0, losses due 
to rounding in the folding process, etc. 

The sensitivity of the targeted search was tested by the detec- 
tion of known faint pulsars in globular clusters, in the framework 
of the targeted search in globular clusters (see Sect. 15.21 ). TableQ] 
gives the results for five pulsars, each observed with three or four 
successive integrations of 160 s (see more details in Sect. 14.21 ). 
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Fig. 1. Signal-to-noise ratio for single observations of PSR 
B1937+21 and reduction in "survey mode". Top: 80-s pointings 
with second setup (24-MHz bandwidth). There were no obser- 
vations in the second semester of 2005 because of a backend 
malfunction. Bottom: 60-s pointings with third setup (48-MHz 
bandwidth). These regular observations of the pulsar were made 
to check the functioning of the backend and monitor the sensi- 
tivity of the setup. 



4.1. Performance validation for "survey mode" 

We regularly observed the millisecond pulsar PSR B1937+21 
(P=l. 557 ms, DM = 7 lcirT 3 pcj3 in order to check the operation 
of the system and to obtain estimates of the achieved detection 
sensitivity. With the last setup, we obtained 99 60-s observations 
on 32 different days (see Fig. [TJ. The average detection SNR is 
28.3 and the median value is 27. 

Our results imply that PSR B 1937 +21 would be detected 
with an SNR of 46 (or 44 when using the median value) with 
an integration of 160 s. Assuming that the main pulse of PSR 
B1937+21 is a source with a mean flux density of 10 mJy, we 
infer a 7-sigma detection threshold of 1.52 mJy (or 1.59 mJy 
when using the median). Taking a safety margin, we can state 
that the sensitivity of the blind survey is about 1.7 mJy for one 
160-s integration with the backend installed in July 2006. 

4.2. Performance validation for "cluster mode" 

We have validated the strategy for processing consecutive obser- 
vations of an accelerated pulsar by applying it to the binary pul- 
sar PSR B 1744-24 A in the globular cluster Terzan 5. This pulsar 
is strong enough (when not in eclipse) to be easily detected and is 
in a close enough binary system to display relatively large varia- 
tions in the period and period derivative within typical observing 
times (a few 160-s observations). We carried out 18 consecutive 
160-s observations on 16 March 2007. 



2 Most data on known pulsars were taken from the ATNF pul- 
sar catalog (http://www.atnf.csiro.au/research/pulsar/psrcat/); see also 
Manchester et al. d2005t 
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Table 1. Known pulsars observed for checking the sensitivity 



Pulsar 


P 


DM 


flux 


N 


obsSNR 


obs flux 




ms 


pc cm" 3 


mJy 






mJy 


PSRB 1620-26 


11.076 


63 


1.6 


4 


19 


1.9-2.1 


PSRB1744-24A 


1.396 


235 


2.5 


3 


30-35 


3.0-3.5 


PSRJ1807-2459A 


3.059 


134 


1.1 


3 


8.4 


0.8-0.9 


PSRB1821-24 


3.054 


120 


1.1 


4 


8 


0.8-0.9 


PSRJ205 1-0827 


4.509 


21 


2.8 


4 


9 


0.9- 1.0 



22 - 
20 - 



survey mode 
1*160sec 



0.8 0.9 1 

Binary phase 
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Fig. 2. Signal-to-noise ratio for observations of PSR B 1744-24 A 
with third setup. Top: single 160-s integrations and reduction in 
"survey mode". Bottom: three 160-s integrations and reduction 
in "cluster mode". 



The 18 observations were then processed in 16 groups of 
three observations, following the procedure described in Sect. 
13.21 The detection SNR values obtained in these analyses are 
shown in Fig. [2] They are perfectly consistent with the coherent 
addition of the SNR values obtained by treating all integration 
separately, except for one point at binary phase 0.98 for which 
the assumption that the pulsar period varies linearly does not 
hold sufficiently over three consecutive observations. 

In our search for new pulsars, we also observed globular 
clusters with known pulsars, which allowed us to check the pro- 
cessing in "cluster mode" (the same validation can also be per- 
formed on field pulsars, but at the cost of some extra observing 
time). Table [TJpresents a summary of the results that we now ex- 
plain in some more details. As an example, PSR J1807-2459A in 
the cluster NGC 6544 has an expected flux density of 1 . 1 mJy at 
1400 MHz, thus would yield an SNR that is quite insufficient for 
detection using only one 160-s integration. Processing of three 
consecutive observations yields an SNR of 8.4, allowing detec- 
tion in a wide search. Similar detections were obtained for sev- 
eral other pulsars. For PSR B 182 1-24 in NGC 6626, t he flux 
density is 0.18 mJy in ATNF, but 1.1 mJy in Lyne et al. d!987l >. 
We obtained an SNR of 8 for four successive 160-s observations, 
which seems to favor the second value. For PSR B 1620-26 in 
NGC 6121, with an expected flux density of 1.6 mJy, we ob- 
tained an SNR of 19 for four 160-s observations. For the field 



Fig. 3. Dependence of the sensitivity on dispersion measure. 
The sensitivity decreases progressively toward higher dispersion 
measures. The diagonal lines correspond to a dispersion repre- 
senting 0.1, 1, and 2 windows per frequency channel, respec- 
tively. The individual points are known millisecond pulsars with 
different symbols depending on their flux density at 1400 MHz. 



eclipsing pulsar PSR J205 1-0827, the expected flux density is 
quite variable and is reported at 2.8 mJy in ATNF, but only 1.1 
mJy in Stappers et al. ([1996 ). We obtained an SNR of 9 for four 
160-s observations, which is consistent with the second value. 
Overall, we can expect to obtain a 8-sigma detection threshold 
at about 1 mJy or slightly below with four consecutive observa- 
tions. 

4.3. Dependence of the actual sensitivity on dispersion 
measure 

The sensitivities quoted above are nominal sensitivities, and the 
actual ones will depend on the dispersion measure at the distance 
of the pulsar. This is true for any processing method. 

In our method, the power spectrum is computed over time 
windows of a small fraction of the tested period and typically 
uses 128 frequency channels. Because of the 1-bit sampling the 
spectrum is normalized and we rely on the variation in the power 
(with time across the windows and with frequency across the 
channels) due to the dispersion of the pulse to detect a pulsar. 
The effect of the dispersion measure on the sensitivity is shown 
in Fig. [3] Starting from the bottom, the first diagonal line is the 
limit below which we can hardly detect pulsars because the dis- 
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Table 2. List of pulsar candidates with SNR> 6.5. 




i cm cmcai)CEinaa> 



canmmnm KcHanimmmD 
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Galactic latitude 

Fig. 4. Map of the observed points in the blind survey in the 
Galactic plane. The green, red, and black open circles corre- 
spond to observations with the first, second, and third setups, 
respectively. The size of the symbols is constant and not repre- 
sentative of the radiotelescope beam, which is elongated in dec- 
lination and increases with declination (thus roughly with lon- 
gitude). The inner rectangle encloses the region of the nominal 
survey. Regions below / = 70° are poorly covered because of 
oversubscription of the radiotelescope time at lower longitudes. 



persion is too low (less than 0. 1 window per frequency channel), 
and the signal is spread equally in all frequency channels. As 
the dispersion increases, the signal for a given window in time 
is spread over fewer frequency channels, and contrast appears 
between channels with and without signal. The sensitivity is op- 
timal (i.e. 0.9 mJy for four integrations of 160 s and reduction in 
"cluster mode") in the region between the first and second diag- 
onal line, i.e., between 0. 1 and 1 window per frequency channel. 
Above the second diagonal line, the sensitivity decreases by a 
factor 2 because the dispersion is between one and two windows 
per frequency channel. Above the top diagonal line, the sensi- 
tivity loss is such that we simply did not search that part of the 
parameter space during processing. These are not strict limits, 
however, at neither the low nor the high end. For example, the 
signal of the pulsar is spread over 20 channels out of 128 for a 
dispersion of 0.05 window per frequency channel, so there is still 
some contrast in the power spectrum between a window with no 
signal and a window with a signal in some part of the band. For 
high dispersion values, the pulse is spread over several windows 
so the fraction of the signal in each window is reduced, and fur- 
thermore the contrast in the power spectrum of a given window 
is also reduced, because the signal will appear in different chan- 
nels for different pulses. 

There are other practical reasons for limiting the search to 
the interval between 0.1 and 2 windows per frequency channel. 
Below 0.1, pulsars with such low dispersion must be nearby and 
should already have been detected, and above 2, pulsars with 
such high dispersion measure are more likely to be distant, thus 
too faint given the decreased sensitivity. 



Name 


RA(B1950) 


Dec(B1950) 


P 

ms 


N 


4091907 


20 00 25.36 


+25 10 


0.876 


10 


2541105 


20 10 18.22 


+37 16 


0.885 


4 


2271416 


20 18 00.63 


+37 16 


1.071 


4 


4830818 


20 29 16.4 


+43 52 


1.047 


7 


3690234 


20 29 22.06 


+38 22 


1.171 


11 


5191628 


20 30 44.2 


+40 56 


0.723 


3 


5141841 


20 33 54.2 


+40 56 


4.347 


3 


4481402 


20 35 56.31 


+41 40 


0.892 


4 


3760145 


20 41 19.79 


+38 44 


1.742 


10 


4942356 


20 41 51.00 


+43 52 


0.997 


5 


4780947 


20 48 30.47 


+43 52 


3.579 


8 


4561319 


20 54 08.7 


+41 40 


2.628 


4 








0.920 




4721049 


21 13 53.4 


+43 08 


2.009 


4 


3630610 


21 36 02.95 


+51 56 


0.954 


4 


3680335 


21 41 02.86 


+54 52 


1.738 


4 


3161822 


22 02 25.0 


+58 54 


2.747 


4 


2791002 


22 07 40.95 


+55 36 


1.720 


4 


3072105 


22 08 50.2 


+58 32 


1.783 


4 


4251900 


22 10 55.6 


+54 08 


0.675 


10 


4012312 


22 16 51.3 


+54 08 


0.952 


7 


3092112 


22 32 51.15 


+58 32 


1.425 


4 


4032330 


22 42 30.55 


+61 28 


1.248 


7 


2210226 


22 49 20.78 


+60 22 


2.108 


4 


2140142 


20 15 12.55 


+36 32 


2.489 


a 


3451149 


22 09 12.42 


+54 30 


1.492 


a 



Notes. <n) not reobserved for lack of observing time 



5. The searches 

Our initial goal was to do a blind survey in the Galactic plane. 
But after two years of unsuccessful searching, we estimated that 
the probability of return per hour of observing time was probably 
higher in targeted searches than in a blind survey, so we decided 
to try them as well. 

There are several possible types of targeted searches for 
millisecond pulsars. The most successful by far is to point at 
Galactic globular clusters. This is where most recent searches 
for such objects have been done (e.g. Possenti et al 2003 
Ransom et al 2005 , Freire et al]2005). Another option is to se- 
lect radiosources with no optical counterpart in catalogs of ra- 
diosources, as done by several groups (Crawford et al. 2000] Han 
et al. 120041 1. and with success by Navarro et al (1995). Pulsars 
can be distinguished from other radiosources in blind surveys by 
their steep spectral index and by a high degree of polarization, 
usually above 5%. Millisecond pulsars tend to show the same, 
sometimes high, degree of linear polarization as ordinary pul- 
sars (Yan et al. 1201 11 1. 



5.1. The blind survey 

The blind survey was conducted in the Galactic plane, where 
most pulsars are likely to be: latitude (b) smaller than 3°, longi- 
tude (/) between 20° and 110°. One could argue that in setting a 
limit in b rather than in height above (or below) the Galactic 
plane, we miss nearby pulsars that are at high b and in the 
Galactic plane. However, as shown below (Sect. 16.21 ). the gain 
in detectable pulsars would be largely offset by the cost in tele- 
scope time. 

A typical observing run was done at a set declination and 
lasted 20 to 120 minutes (depending on the declination and the 
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Table 3. List of the observed globular clusters 



Cluster 


DM 


Np 


date 


N 


Cluster 


DM 


Np 


date 


N 




pc cm -3 




yyyymmdd 






pc cm~ 3 




yyyymmdd 




NGC 5986 


169-171 





20070213 


4 








20070815 


18 


NGC 6093 


108-109 





20070312 


4 








20070816 


18 


NGC 6121 


63 


20 


20070213 


4 


NGC 6553 


237-350 


4 


20070313 


4 


NGC 6171 


83-86 


4 


20061119 


3 








20071211 


8 








20070213 


4 


NGC 6624 


87 





20070427 


4 








20070426 


4 


NGC 6626 


121 


5 


20070211 


4 








20071005 


10 


NGC 6712 


235-333 


3 


20061114 


3 








20071201 


10 








20070313 


4 


NGC 6218 


70-75 


7 


20061119 


3 








20071221 


18 








20070213 


4 


NGC 6717 


163-193 


2 


20061114 


3 


NGC 6254 


73-81 


8 


20061119 


3 








20070923 


18 








20070312 


4 


NGC 6760 


208-310 


1 


20061203 


4 








20070426 


4 


NGC 6809 


103-113 


6 


20061203 


4 








20070624 


18 








20070117 


4 








20070718 


18 








20070504 


18 








20071101 


10 








20071121 


10 


NGC 6256 


418-554 





20071201 


10 


2MASSGC01 


144-292 


13 


20070210 


4 


NGC 6266 


115 


3 


20070313 


4 








20070719 


18 


NGC 6304 


224-312 


4 


20070119 


4 


2MASSGC02 


127-175 


7 


20070211 


4 








20070210 


4 








20070922 


10 








20070815 


18 


FSR1767 


61-110 


25 


20070922 


8 


NGC 6342 


103 





20070426 


4 


Glimpse-C02 


294-560 


5 


20080919 


4 


NGC 6366 


69-85 


13 


20070119 


4 








20081122 


5 








20070516 


4 


IC 1276 


176-251 


5 


20070427 


4 


NGC 6440 


224 





20070427 


4 








20071113 


18 


NGC 6441 


366-451 





20070427 


4 


TerzanlO 


268-413 


5 


20071211 


8 


NGC 6544 


134 


16 


20061114 


3 


Terzanl2 


197-332 


7 


20061114 


3 


NGC 6540 


212-315 


13 


20070119 


4 








20070211 


4 








20070516 


18 








20071012 


18 



alloted time), shifting the telescope (in fact the receiver, since 
this is a transit instrument) by one beam in right ascension after 
each short integration. A better observing strategy was to shift by 
multiples of one beam and return to that declination at later dates 
to fill the gaps. The selected declinations were also separated by 
one beam, so that all pointings were separated by at least one 
beam in either direction. 

The total amount of telescope time attributed to this project 
over the years was 549.3 hours, at which point the time alloca- 
tion committee wisely decided to cut its losses. We observed 186 
points in the sky with the first setup (sensitivity ^3.5 mJy for 
integrations of 240 s), 2829 points with the second setup (sensi- 
tivity =i 2.2 mJy for integrations of 200 s) and 2542 points with 
the last setup (sensitivity ^1.7 mJy for integrations of 160 s). 
Seventy points were observed several times, leading to a total of 
5487 different points observed out of a total of 22057, that is, 
a coverage of 25% of the planned survey. The actual sensitivity 
was different in some of the pointings. In the second semester 
of 2004, the sensitivity progressively degraded from 2.2 to 2.7 
mJy due to an unnoticed desynchronization in the acquisition. 
A bug in the treatment software in late December 2004 low- 
ered the sensitivity by about half for 57 points. We attributed 
a slightly higher sensitivity (by an arbitrary amount of 20%) to 
the 70 points observed repeatedly, because of the better chance 
of detecting pulses as the effective flux may vary and as the fixed 
set of tested values in the parameter space may better match the 
actual parameters of the pulsar. 

A map of the survey region in Galactic coordinates with the 
covered region shown as circles (colors depending on setup) is 
presented in Fig. |4] The figure shows that a few points were ob- 
served outside the region of the survey, either by mistake, or by 
opportunity (to save repointing the telescope). The figure also 



shows that regions below I = 70° are poorly sampled, because 
of oversubscription of the telescope time by a factor greater 
than 2 at the corresponding right ascensions. While millisecond 
pulsars are likely to be more frequent at lower Galactic longi- 
tudes (closer to the Galactic center), their a priori detectability 
is strongly reduced by pulse smearing at longitudes below 50° 
(see Sect. |6]and Fig. |5}- However, the drawback of higher decli- 
nations (and thus higher longitudes) is that the antenna aperture 
efficiency is progressively reduced, with a loss of sensitivity of 
15% at I a 85° and 28% at I * 1 10° (Fouque et al. [19901 About 
half of our pointings suffer a loss of sensitivity of 20% or more 
because of that effect. 

During the survey we identified 25 pulsar candidates with an 
SNR^7. They are listed in Table [2] together with their coordi- 
nates (1950 equinox), probable period (in ms) and the number 
of 160-s integrations in the reobservations. Note that the last two 
candidates were not reobserved for lack of observing time. Their 
reobservation and analysis in "cluster mode" (see Sect. 13.2b did 
not enable us to confirm them. Some of these presumably false 
detections can be explained by interferences from radars at civil- 
ian and military airports near Paris. 

5.2. Globular clusters 

The central regions of globular clusters are favorite targets when 
searching for millisecond pulsars, because they are regions of 
very high stellar density. In such regions, multiple encounters 
between stars favor the formation of binary stars (Baiivn [T9951 l. 
and, when one of the stars is a pulsar, favors the acceleration of 
the pulsar to periods of about a millisecond (pulsar recycling). 
Furthermore, a large number of binary stars in globular clusters 
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may in fact be primordial binaries (Hut et al. 1992), thus in- 
creasing the probability of finding recycled pulsars in globular 
clusters. 

We selected globular clusters as targets using three criteria: 

- Avoid the zone accessible to the radiotelescope in Arecibo, 
which, with its large collecting area and its new ALFA re- 
ceiver, can conduct much more sensitive searches than us. 

- Aim at nearby globular clusters, with high stellar density, 
and where no millisecond pulsars have been detected so far. 

- Select preferentially those clusters where the dispersion 
measure is highest, and hence where pulsars are difficult to 
detect with some of the search techniques used in past sur- 
veys (see below). 

To test the feasability of our project, we selected in the ATNF 
the 32 pulsars discovered in globular clusters as of April 2006 
with a published flux density. If the flux density was given at 
400 MHz, it was converted to 1400 MHz assuming a spectral 
index of -1.8. We then calculated the number Np of these 32 
pulsars that would be detected if they were all in each of the 
different Galactic globular clusters. We selected those clusters 
in which we have a reasonable chance of detected millisecond 
pulsars taking Np, the sensitivity of our search, and the value of 
the dispersion measure at the distance of the cluster into account. 
Other clusters were added in the course of the search to test the 
sensitivity for reduction in "cluster mode" (see Sect. 14.21 ) or to 
optimize the use of the available telescope time. 

We obtained 28.33 hours (+ 3 of discretionary director's 
time) of telescope time for this project. The list of observed glob- 
ular clusters and their characteristics are given in Table [3] The 
successive columns are the name of the cluster, the range in dis- 
persion measure (from the model of Cordes & Lazio 2002]), or 
the dispersion measure (from the ATNF catalog) if millisecond 
pulsars have been detected in the cluster, in cirT 3 pc, the expected 
number Np of pulsars with a flux density above 1 mJy (assuming 
that all 32 know millisecond pulsars are in that cluster), the date 
of observation, and the number of pointings. 

Candidates were detected in the clusters NGC 6171, NGC 
6254 and NGC 6809, which were reobserved repeatedly. No new 
millisecond pulsar was detected. 

5.3. Bright polarized point-like radiosources with steep 
spectra 

In the second targeted search, the candidate pulsars were se- 
lected from the catalog of steep spectral index radiosources pub- 
lished by De Breuck et al. (f2p00). While their catalog was con- 
structed with the aim of obtaining a large sample of high-redshift 
radiosources, because spectral index tends to increase with the 
redshift of optically identified sources, a large negative spectral 
index is also a property of millisecond pulsars. 

We extracted the information on polarization and on source 
size from the NVSS0 catalog of radiosources. From this com- 
plete list, we extracted the sources that were unresolved and po- 
larized (> 5%). These criteria give 26 sources, two of which are 
in fact known pulsars (PSR J0630-2834 and PSR J2313+4253). 
We removed them from the list and checked that the 24 other 
sources are not in any database (Simbad, NED, Hyperleda). The 
final list is given in Table [4] together with their 1400-MHz flux 
density in mJy, spectral index between 1400 and 326 MHz, per- 
cent polarization, the date of our observation, and the number 

3 http://www.cv.nrao.edu/nvss/ 



Table 4. List of observed bright polarized unidentified ra- 
diosources with steep spectrum 



Name 


S1400 


a 


% pol. 


date 


N 




mJy 






yyyymmdd 




WNJ075 1+3300 


10.7 


-1.47 


5 


20080107 


18 


WN J0829+3834 


10.9 


-1.49 


5 


20080202 


18 


TN J0856-1510 


68.5 


-1.41 


12 


20080304 


18 


WN J0928+6003 


10.4 


-1.63 


7 


20080402 


18 


WN Jl 101+3520 


14.9 


-1.32 


5 


20080204 


18 


WN Jllll+3311 


12.3 


-1.32 


6 


20080319 


18 


WN Jl 152+3732 


14.7 


-2.18 


5 


20080110 


18 


WN J1258+3212 


18.3 


-1.52 


8 


20080306 


18 


WN J 1300+53 11 


18.1 


-1.44 


5 


20080214 


18 


WN J1330+6505 


9.5 


-1.34 


5 


20080413 


3 


WN J1337+3401 


13.2 


-1.39 


7 


20080116 


12 


WN J 1359+7446 


11.2 


-1.92 


5 


20080413 


4 


WN J 1400+4348 


18.2 


-1.49 


6 


20080205 


4 


WN J1421+3103 


11.6 


-1.43 


5 


20080413 


3 


WN J1435+3523 


20.8 


-1.35 


7 


20080320 


12 


TN J1515-2651 


56.7 


-1.34 


9 


20080218 


18 


WN J1606+4142 


10.9 


-1.31 


10 


20080407 


18 


WN J1723+5822 


10.1 


-1.46 


11 


20080111 


8 


WN J1734+3606 


7.1 


-1.31 


14 


20080117 


8 


WN J1739+5309 


9.8 


-1.32 


5 


20080111 


8 


WN J1819+6213 


11.6 


-1.45 


7 


20080123 


12 


WN J1923+6047 


11.8 


-1.30 


8 


20080419 


18 


TN J2028-1934 


131.4 


-1.34 


7 


20080302 


15 


WN J2221+3800 


13.2 


-1.32 


5 


20080422 


18 



of 160-s integrations. More information on the sources can be 
found in De Breuck et al. (120001 ) . 

We obtained 18.67 hours of telescope time for this project 
and observed all the sources in the list. None of them turned out 
to be a pulsar. 

5.4. Faint polarized point-like radiosources 

In the third search, we extracted the candidate pulsars from the 
NVSS catalog. This catalog contained 97 of the 1775 pulsars in 
the ATNF catalog at the time of the project (early 2008). From 
the NVSS catalog, we first extracted the sources with 1400-MHz 
flux densities between 2 and 10 mJy. We did not try to determine 
their spectral index, as in the previous set of targets, because 
there is no hope of finding these faint sources in catalogs at other 
frequencies. 

We then selected those sources that were polarized (po- 
larized flux density > 20% of total flux density) and with 
SNR in total and polarized flux density above 3. This pro- 
duced a list of 2105 sources. In this list, we found six pul- 
sars (PSR J0358+5413, PSR J0538+2817, PSR J0814+7429, 
PSR J0922+0638, PSR J1745-3040, PSR J1840+5640), none 
of which has a very short period. But one should note that 
none of the 45 known millisecond pulsars with periods less than 
3.3 ms has 1400-MHz flux densities above 2 mJy, except PSR 
B1937+21, so they cannot be in that catalog. We then further 
restricted the sample to the Galactic plane {\b\ < 12°), also 
eliminating the region 75° < / < 165° and \b\ < 5.5° (re- 
cently surveyed by Hessels et al. 2007), the region of our own 
blind survey, and the right ascensions heavily oversubscribed 
at Nangay (17h00-19h40). We removed all sources identified in 
NED, SIMBAD, and 2MASX. We also carefully examined the 
postage-stamp images of all the sources and removed those that 
looked like noise, as well as close double sources (which are 
likely to be double-lobed radiogalaxies). This brought the to- 
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tal number of targets to 92. Some ambiguous cases (multiple 
sources, not quite point-like sources) were left in the sample, 
since at this low flux level it is sometimes difficult to distinguish 
real features from noise. 

We were given 2 1 .67 hours of telescope time for this project, 
which enabled us to observe 65 of them (with one source ob- 
served twice). They are listed in Table[6j which gives the NVSS 
name, the flux density at 1400 MHz in mJy, and the percent po- 
larization. This last targeted search did not produce any new mil- 
lisecond pulsar, which seems to confirm that most radiosources 
with ultrahigh polarization are radio-loud active galactic nuclei 
(Shietal.|MS. 

The conclusion of the two targeted searches in catalogs of ra- 
diosources is that the luminosity function of millisecond pulsars 
is simply too faint to allow a meaningful number of them to be 
found in large surveys like NVSS or FIRST. The small number 
of sources observed does not allow us to be more quantitative. 

6. The Galactic population of millisecond pulsars 

We now examine whether the null result of our blind survey 
is statistically compatible with the Galactic distribution of mil- 
lisecond pulsars as we know it. The method is the following. 
We adopt a model for the distribution of millisecond pulsars in 
space, rotation period, and luminosity, produce a large number 
of Galactic populations of synthetic millisecond pulsars, and, for 
each realization of the simulation, we compute how many syn- 
thetic pulsars should be detected in our blind survey, by com- 
paring the properties of each synthetic pulsar (position, rotation 
period, flux density) with the coordinates and actual sensitivity 
of each of our 5487 pointings, taking the effect of the disper- 
sion measure on the detectability into account in each case. The 
number of detected pulsars should be zero in most realizations, 
otherwise our model overestimates the true number of millisec- 
ond pulsars. 

6.1. Model for the Galactic population of pulsars 

Several models have been proposed over the years to charac- 
terize the population of pulsars in the Galaxy, either ordinary 
(Lyne et al. I1985I Yusifov an d Kiigi ik, I2004I Lorimer et al. 
2006 ; Faucher-Giguere & Kaspi 2006), millisecond (Lorimer et 
al. [19961 C ordes & Chernoff, [19971 Story et al. I20071 . or both 
(Lyne et al. [19981 Faucher-Giguere & Loeb,HJTID- 

The population of millisecond pulsars is too small to allow 
for a precise modeling of their spatial distribution. However, they 
are thought to have larger radial and vertical scale lengths than 
those of ordinary pulsars, because they are likely to have lost the 
memory of their birth location due to a high initial velocity. 

We adopt the spatial distribution used by Faucher-Giguere 
& Loeb (2010). The radial density in the Galactic plane is 
p(r) oc exp(-r 2 /2cr 2 ), where r is the radial distance from the 
Galactic center and cr = 5 kpc. For ordinary pulsars, cr — 3.6 kpc 
(Yusifov and Kiigiik, 2004). Our simulations and others (Lorimer 
et al. 1 19961 1 have shown that the exact radial dependence is in 
fact not important, because most millisecond pulsars are faint, 
and the heavily populated region toward the Galactic center is 
"obscured" by a high dispersion measure (see Fig.|5]l. Most de- 
tectable millisecond pulsars are thus within 2 or 3 kpc of the 
Sun (which in this paper is assumed to be located 8 kpc from the 
Galactic center). The distribution perpendicular to the Galactic 
plane is a simple exponential : p(z) oc exp(-\z\/h), where z is the 
distance from the Galactic plane and the scale height h = 500 



pc (Cordes & Chernoff [T9971 Story et al. EOOTY This value is at 
the upper end of those estimated for ordinary pulsars (330-600 
pc). Adopting a higher value for the scale height slightly reduces 
the number of pulsars at low Galactic latitudes, and hence the 
probability of detection in our blind survey. We do not take into 
account a possible increase in h with r (Yusifov 2004), which 
would reduce the number of pulsars in the Galactic plane even 
more. 

The total number of millisecond pulsars is highly uncertain. 
Lorimer et al. (1995) estimate that there are 40 000 millisecond 
and low-mass binary pulsars above L„„„ = 2.5 mJykpc 2 at 436 
MHz. Lyne et al. (1998) estimate that there are 30 000 millisec- 
ond pulsars above L mi „ = 1 mJykpc 2 at 400 MHz. Desvignes 
(20092009) finds 46 000 millisecond pulsars above L mi n = 1 
mJykpc 2 at 1400MHz using model C of Lorimer et al. (2006), 
which amounts to 23 000 above 0.2 mJykpc 2 with our adopted 
luminosity function. We adopted a total number of 30 000 mil- 
lisecond pulsars above L ml „ = 0.2 mJykpc 2 at 1400 MHz, which 
is a happy medium between the above estimates, using a spectral 
slope of -1.8 to convert L mi „ to the right frequency. 

The observed distribution of millisecond pulsar rotation peri- 
ods P is a power law of index +3 (Hessels et al. 2007). However, 
the paucity of fast-spinning pulsars is the result of selection ef- 
fects (see Sect. Q] and Lorimer et al. [1996). Following Lorimer 
et al. (1996), we adopt a distribution of rotation periods that is 
uniform in log(f). The lower limit of 0.65ms is the lowest theo- 
retically possible one (Haensel et al. 1995), and the upper limit 
reflects an observed cutoff. Cordes & Chernoff (1997) assume a 
period distribution in the form P l . For the rotation periods of 
ordinary pulsars, a log-normal distribution is generally assumed 
(Lyne et al. [19831 Lorimer et al. 120051 . 

The distribution of luminosities is also modeled by a power 
law, down to the lowest luminosity (Lorimer et al. 120071 Hessels 
et al. [20071 Bagchi & Lorimer [20T0JI: p(L) = L" if L > L mi „, 
and otherwise, with n around unity. We adopted n = 1. The 
luminosity distribution of ordinary pulsars has been modeled by 
one or two power laws (Lyne et al. 1 19981 Guseinov et al. 2003). 
Models where the luminosity distribution depends on the pulsar 
period and its derivative have also been considered, for both mil- 
lisecond and ordinary pulsars (Story et al. 2007). The flux den- 
sity was computed by multiplying the luminosity by the square 
of the distance to the synthetic pulsar. 

Finally, we take possible pulse smearing by the interstellar 
medium into account, using the Galactic free electron model 
NE2OO10 of Cordes & Lazio d20021 l to estimate the disper- 
sion measure (DM) at the distance of the synthetic pulsars. If 
DM > 86.38P, where P is the period in ms and DM is in cm~ 3 pc, 
the sensitivity is reduced by half. If DM > 172.76P, the sensi- 
tivity is reduced by a factor 3 or more, and we consider that the 
pulsar is undetectable (see Fig. [3j. On the other hand, we ignore 
the effects on the flux density of diffractive interstellar scintilla- 
tion, which can make the measured flux density vary by a large 
factor, as shown for PSR B1937+21 in Fig.Q] 

6.2. Results of the simulations 

We did many simulations of the distribution and detectability of 
the pulsar population in the region of our blind survey, using the 
above model and parameters. Out of a population of 30 000 syn- 
thetic pulsars in the Galaxy, only about 4600 short-period mil- 
lisecond pulsars (P < 3ms) fall in the region of the survey ex- 
tended to 0° < I < 120° and \b\ < 4° for comparison with Fig. [4] 

4 http://rsd-www.nrl.navy.mil/72 1 3/lazio/ne_model 
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Fig. 5. One simulation of the population of 30 000 millisec- 
ond pulsars in the Galaxy. Top: Distribution of all (4589) short- 
period pulsars in the region of the blind survey and correspond- 
ing histogram of flux densities; the solid curve is the mean his- 
togram from 1000 simulations. Bottom: distribution of the (730) 
pulsars remaining after eliminating the undetectable ones (be- 
cause of pulse smearing at high dispersion measure) and corre- 
sponding histogram of flux densities. Both histograms are lim- 
ited to flux densities higher than 0. 1 mJy for clarity. 



There are about 5700 long-period millisecond pulsars (P > 3ms) 
in that region, which thus contains a bit more than a third of the 
total population of synthetic pulsars. We show the result of one 
simulation in Fig. [5] The distribution of the short-period pul- 
sars is shown in the top panel, together with the histogram of 
flux densities above 0. 1 mJy. The smooth curve is the mean his- 
togram for a thousand simulations. The pulsars are concentrated 
in the region nearest to the Galactic center. 

However, when considering the effect of high DM, that is, 
eliminating all the pulsars above the top diagonal line of Fig.[3j 
the overall population of detectable pulsars in the region of the 
survey is reduced by a factor 6. This is shown in the bottom panel 
of Fig. [5] where the region near the Galactic center is no longer 
heavily populated, and the region 40° < I < 80° becomes the 
most populated. This depopulation of the detectable pulsars is 
not so critical if one only considers the pulsars with a high flux 
density, say above 1 mJy. This population is reduced from 17 to 
9, that is, by about half, because the depopulation affects mainly 
pulsars with a high DM, which are necessarily more distant, thus 
fainter and more difficult to detect. 

The number of detectable pulsars for a given sensitivity can- 
not be read directly from the histogram of the lower panel of 
Fig. [5] because the actual sensitivity will depend on the DM at 
the distance of the pulsar. Taking the example of our project, 
we find that 4.2 pulsars in the extended region of our survey 
are detectable for the nominal sensitivity of 1 .7 mJy of our last 
setup (taking the effect of large DM into account), and 7.7 pul- 
sars are detectable for a sensitivity of 0.9 mJy (which is achieved 
with four integrations of 160 s and reduction in "cluster mode"). 
This means that, if we want to double our chances of detecting a 



pulsar, we have to increase our sensitivity twofold. This can be 
achieved by a larger bandwidth, but, if keeping the equipment 
unchanged, it would require increasing the telescope time four- 
fold, which is not an economical strategy. It is more efficient to 
extend the region of the survey rather than increase the sensi- 
tivity (provided the added region is as rich in potential pulsars 
as the initial region). Only at much higher sensitivities does the 
number of synthetic pulsars go up significantly with increased 
sensitivity (see histograms of Fig. [5). 

Now, to estimate whether a synthetic pulsar is detected in any 
one of our pointings, we model the beam pattern of the radiote- 
lescope by a sine 2 function in both coordinates (modeling by a 
Gaussian is only acceptable within the half -power beamwidth), 
and compute the sensitivity correction at that pointing as As = 
sinc 2 {da)sinc 2 {dd), where da and d5 are the angular distances 
of the synthetic pulsar to the center of the beam in right ascen- 
sion and declination, respectively. The half -power beamwidth is 
3.6'x22' at 1400 MHz and zero declination, and there is an ad- 
ditional correction in declination above 22.7°. The beamwidth is 
about a factor 1.32 larger in declination at our highest declina- 
tion. We also consider the loss of gain due to variable aperture 
efficiency and atmospheric absorption, using the equation given 
by Fouque et al. (1990, their equation 1). 

The results of 200 000 simulations are given in Table [5] 
where N p is the number of detected synthetic pulsars, P s the 
probability of detecting N p pulsars in our survey of 5487 point- 
ings, Pi 7 the probability of detecting N p pulsars if we had com- 
pleted the whole survey with a nominal sensitivity of 1.7 mJy 
from beginning to end, and Po.9 the probability of detecting N p 
pulsars if we had surveyed the more densely populated region 
\b\ < 2° and 45.0° < I < 68.0° (3761 pointings) with a higher 
sensitivity of 0.9 mJy (by doing four 160-s integrations at each 
pointing and reducing the data in "cluster mode"). 

In summary, we had one chance in four of detecting one or 
more short-period pulsars in our incomplete survey. Were we 
to redo the complete survey with the sensitivity reached in the 
last setup, we would have a 79% chance of success, but even 
so the harvest would be far from spectacular. And this ignores 
the effect of interstellar scintillation, which reduces our chances 
of success by an unknown amount. While the third option of 
surveying a smaller region with a higher sensitivity may a priori 
seem attractive, it has a low return and is in practice unfeasible 
at a transit instrument like the Nancay radiotelescope, because 
it would monopolize the right-ascension slot 19h-20h for about 
two years. 

As mentioned in Sect. 15. II we miss nearby high-latitude mil- 
lisecond pulsars when limiting the blind survey to low latitudes. 
If we set no limit in b (there is of course no way of limiting the 
survey in linear height above the Galactic plane), we find that 
there are about 6510 (rather than 4600) short-period millisec- 
ond pulsars in the region of our survey, 1980 (rather than 730) 
detectable ones, and 29 (rather than 9) above 1 mJy. However, 
we would have to scan about half the sky visible at Nangay to 
complete the survey, so this is another unrealistic option. 

Comparing our simulations with the population of known 
millisecond pulsars, we find 14 known field millisecond pulsars 
with a known flux density at 1400 MHz in the region of our blind 
survey, and the faintest one has a flux density of 0.05 mJy. Down 
to that value, there are about 33 synthetic millisecond pulsars, so 
that about half the total population predicted by our simulations 
in that region have already been detected. 

Incidentally, the sensitivity of our blind survey is such 
that, among these 14 known millisecond pulsars, only PSR 
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Table 5. Probability of detecting short-period pulsars 





Ps(%) 









0.763 


0.213 


0.697 


1 


0.151 


0.223 


0.174 


2 


0.055 


0.184 


0.073 


3 


0.020 


0.145 


0.035 


4 


0.006 


0.099 


0.014 


5 


0.003 


0.059 


0.004 


6 


0.0006 


0.038 


0.0014 


7 


0.0002 


0.020 


0.0006 



B1937+21 and PSR J2007+2722 (which has a period of 24 ms) 
can be detected with our setup. 

To put the results of our blind survey in perspective, we men- 
tion the results of another blind survey for millisecond pulsars 
conducted at Nangay at about the time that we started our own 
survey. A survey of 45000 pointings in the region 180° > / > 
-15° and \b\ < 3° was conducted in 1997-98 (Foster et al. fTWTt , 
using a different backend and method of analysis. The frequency, 
bandpass, integration time, and nominal sensitivity of the survey 
were 1360 MHz, 144 MHz, 120 s, and 0.6 mJy, respectively, 
and the sampling resolution was 1 .5 MHz in frequency and 60 
fis in time. It produced only two new pulsars, both with periods 
longer than 100 ms (Ray et al. |1999l >. The reanalysis of 15% of 
the stored data with a different method produced 341 candidate 
millisecond pulsars, and observing 200 of them again did not 
reveal any new pulsar (Desvignes 2009). With our simulations, 
we find that there are 30 synthetic millisecond pulsars in the re- 
gion of the survey above the flux limit, but only 14 of them are 
detectable with our processing method in the best-case scenario 
(i.e. if the telescope points exactly at the pulsar in each case). 
This confirms what was probably widely believed, that it is no- 
toriously difficult to detect millisecond pulsars. As one referee of 
our telescope time requests stated, it is a high-risk, high-reward 
project. 

7. Conclusions 

The present project was started with no assumptions on the to- 
tal population of millisecond pulsars (except that it must be 
large) and with the specific aim of detecting a population of fast- 
spinning (sub)millisecond pulsars. The results of the small sur- 
veys targeted at unidentified radiosources suggest that, despite 
PSR B1937+21, bright fast-spinning pulsars are rare. The null 
result of the blind survey agrees with our present understanding 
of the field population of millisecond pulsars, the order of mag- 
nitude of the total population, its absence of strong clustering, 
and the distribution of periods. Our conclusion is that the present 
models of the Galactic population of millisecond pulsars are not 
unduly pessimistic and that the population of fast-spinning mil- 
lisecond pulsars cannot be unexpectedly large. 

This project was also fruitful in a broader context: over 24 
students worked as interns on the project over the years, and 
most of them subsequently made use of the acquired skills in 
their career. 

Finally, we hope that future prospective searches for mil- 
lisecond pulsars will successfully make use of our experience 
and of the numbers given here. 
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Table 6. List of observed faint polarized point-like radiosources 



Name 


S1400 


% pol 


date 


N 


Name 


S1400 


% pol 


date 


N 




mjy 




yyyymmdd 






mjy 




yyyymmdd 




NVSS J002145+744454 


9.7 


26.0 


20080905 


3 


NVSS J061712+260152 


7.5 


27.5 


20081103 


4 


NVSS J003706+5 15737 


7.3 


40.7 


20080905 


3 


NVSS J070822-151518 


10.0 


26.1 


20081204 


3 


NVSS J004419+514533 


2.9 


84.5 


20080905 


3 


NVSS J071545+051635 


8.5 


24.9 


20081012 


4 


NVSS J012635+704215 


4.4 


76.8 


20080905 


4 


NVSS J071933-365933 


7.6 


26.4 


20081106 


3 


NVSS J012705+703605 


3.5 


75.1 


20080909 


3 


NVSS J072438-271618 


5.0 


43.0 


20081106 


4 


NVSS J0128 18+704655 


5.3 


64.0 


20080909 


3 


NVSS J07250 1+004534 


7.5 


41.6 


20081204 


4 


NVSS J012932+701545 


4.2 


50.5 


20080905 


4 


NVSS J073949+004434 


7.1 


28.6 


20081106 


4 


NVSS J013201+704729 


3.3 


113.3 


20080909 


4 








20081204 


4 


NVSS J030809+5 12443 


7.4 


32.7 


20081003 


3 


NVSS J082008-363214 


6.3 


34.3 


20080703 


3 


NVSS J034707+402127 


9.4 


28.2 


20081003 


4 


NVSS J082838-380231 


3.4 


68.2 


20080803 


4 


NVSS J035429+643653 


9.5 


22.5 


20081001 


3 


NVSS J083035-3 12745 


6.5 


32.8 


20080703 


3 


NVSS J042204+363943 


9.3 


30.6 


20081001 


3 


NVSS J090355-384439 


3.0 


67.3 


20080703 


4 


NVSS J043739+301858 


5.9 


66.3 


20081001 


4 


NVSS J155955-382320 


6.0 


35.3 


20080911 


3 


NVSS J050224+250015 


5.1 


43.3 


20081001 


4 


NVSS J1601 16-381526 


5.4 


68.1 


20080805 


3 


NVSS J052352+164533 


5.2 


38.7 


20081210 


3 


NVSS J162707-364109 


3.9 


55.1 


20080805 


3 


NVSS J052906+215305 


6.7 


59.7 


20081212 


3 


NVSS J162943-340416 


8.4 


27.3 


20080911 


3 


NVSS J053220+ 163421 


8.8 


22.8 


20081103 


3 


NVSS J163550-375558 


4.7 


54.9 


20080706 


3 


NVSS J053959+193737 


5.4 


37.6 


20081113 


3 


NVSS J164751-341010 


4.2 


48.3 


20080805 


4 


NVSS J053959+232629 


2.5 


112.4 


20081011 


3 


NVSS J165146-353221 


5.1 


39.8 


20080706 


4 


NVSS J054134+302901 


9.5 


22.4 


20081210 


4 


NVSS J165753-305936 


5.1 


52.2 


20080911 


4 


NVSS J05433 1+500635 


4.1 


49.5 


20081212 


4 


NVSS J194518+101520 


9.9 


20.2 


20080809 


3 


NVSS J054456+491607 


5.0 


47.2 


20081103 


4 


NVSS J 1950 16+404408 


4.9 


44.9 


20080807 


3 


NVSS J054532+215524 


8.7 


24.7 


20081011 


3 


NVSS J200 102+420928 


4.6 


100.7 


20080807 


3 


NVSS J054650+495 147 


6.3 


44.3 


20081113 


4 


NVSS J200109+1 12314 


8.4 


26.8 


20080809 


3 


NVSS J054747+392734 


9.3 


23.9 


20081210 


4 


NVSS J202139+233729 


4.1 


54.1 


20080809 


4 


NVSS J055 120+032729 


8.2 


26.3 


20081212 


4 


NVSS J203453+214620 


4.6 


88.7 


20080807 


4 


NVSS J055600+185916 


9.1 


27.5 


20081103 


4 


NVSS J204040+304747 


3.5 


62.9 


20080711 


3 


NVSS J060348+215030 


3.6 


57.2 


20081011 


4 


NVSS J210904+382547 


7.8 


26.0 


20080711 


3 


NVSS J060348+343522 


10.0 


22.0 


20081113 


4 


NVSS J21 1354+61 1530 


3.7 


61.4 


20080711 


4 


NVSS J06041 1+31 1018 


9.7 


21.8 


20081212 


4 


NVSS J215324+395318 


8.6 


24.0 


20080711 


4 


NVSS J060820+171726 


8.0 


26.0 


20081210 


4 


NVSS J223844+45 1415 


10.0 


20.4 


20080708 


3 


NVSS J06101 1+150106 


7.4 


28.8 


20081011 


4 


NVSS J230058+713519 


9.8 


26.4 


20080708 


3 


NVSS J061040+332755 


8.6 


24.1 


20081113 


4 


NVSS J234546+544241 


8.2 


24.9 


20080708 


4 
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